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ABSTRACT 



High-resolution Hubble Space Telescope (HST) archive imaging and high- 
dispersion spectroscopy are used to study the complex morphological and kine- 
matical structure of the planetary nebula, NGC 2818. We analyze narrow band 
Ha, [O III], [N II], [S II] and He II images, addressing important morphologi- 
cal features. Ground-based longslit echelle spectra were obtained crossing NGC 
2818 at five different positions to precisely determin e kinematical features in 



the s tructure of the nebula. A distance of 2.5 kpc (I Van de Steene &: Zijlstra 



19951 ) was used to determine physical scales. Constructing models to fit the 
data with modern computational tools, we find NGC 2818 is composed by: (1) 
a non-uniform bipolar structure with a semi-major axis of 0.92 pc (75arcsec), 
possibly deformed by the stellar wind, (2) a 0.17pc (14arcsec) diameter central 
region, which is potentially the remnant of an equatorial enhancement, and (3) 
a great number of cometary knots. These knots are preferentially located in- 
side a radius of 0.24 pc (20arcsec) around the central star. The major axis of 
the main structure is oriented at « ~ 60° with respect to the line-of-sight and 
at PA = +89° on the plane of the sky. Expansion velocities of this nebula are 
Vpoi = 105kms~ and Veq = 20kms~ , which lead to our estimate of the kine- 
matical age of Tk ^ 8,400 ± 3,400yr (assuming homologous expansion). Our 
observations do not support the idea that high velocity collimated ejections are 
responsible for the formation of microstructures inside the nebula. We determine 



the systemic velocity of NGC 2818 to be I/rel = +26 ± 2kms"\ 

Subject headings: planetary nebulae: individual: NGC 2818 — ISM: kinematics and 
dynamics 
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Introduction 



1.1. Morphology of Planetary Nebulae 



The Planetary Nebula (PN) phase occurs late in the evolution of low-mass stars 
when UV radiation from the star produces a photoionization-recombination balance in an 
expanding gaseous shell which had been expelled during the asympototic giant branch 
(AGB) phase. Supersonic expansion velocities of 30-40 kms~^ are seen in many planetary 
nebulae (PNe) shells, in which the sound velocity is ~ l Okms"^. To explain why th e PNe 
do not follow a free thermal expansion (sound velocity), iKwok. Purton fc Fitzgerald (119781 ) 
proposed the Interacting Stellar Wind (ISW) model, whereby a denser and slower AGB 
wind interacts with a new fainter and faster stellar wind. The existence of a fast wind 



around the central star of PN was confirmed years later b y means of UV o 
the International Ultraviolet Explorer (lUE) satellite (e.g. 



Heap et al. 



Dservations with 



19781). 



Improvements in optical imaging have revealed that an understanding of the 
morpholog y of PNe requ ires more complex explanations and depends on additional physical 



processes. 



Balickl (Il987l ) proposed a modification to the ISW model, the Generalized 



Interacting Stellar Wind (GISW) model, which includes an equatorial density enhancement 
to explain elliptical and bipolar morphologies. This mo del \vas sub seque ntly confirmed as 
an ap propriate approximation and further developed by 



Ickd mm and 



Frank fc Mellema 



(119941 ). Since then, many other morphological features have been dis covered in deeper and 



better resolved observation s: open bipolar lo 



multi polar struct ures (e.g. 



2008 ). rings (e.g. 



Lopez et al. 



Vazquez et al. 



2002 



) es (called 'diabolos'; e.g 



1998 



Guerrero fc Manchado 



Bond et al 



2003 



Monteiro et al 



19981: 



MitcheU et al 



2000), 



Vazquez et al. 



2006 



Jones et al. 



2010l ). and microstructur es (also called Low Ionization Structures or LIS) such as ansa e 



knots, jets and filaments fjGoncalves. Corradi fc Mampaso 



2001 



Miszalski et al. 



2009H ). In 



some cases, point-symmetry appears as an additional ingredient to the general morphology 
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(e.g. iGuerrero. Vazquez &: Lopea Il999t ISabin et al.l 1201 ll ). As PNe have different 
observed morphological features, many PNe cannot clearly be assigned membership to a 
morphological class based on a single defining feature (round, bipolar, or elliptical). 



Bipolarity has been explained since the pioneering work of iBalickl (119871 ) as a 
consequence of an equatorial density enhancement of material formed by a mass-loaded 
'superwind' ejected in the late-AGB phase. However, it is still poorly understood how 
this superwind is trig gered and why i t diverges from a spherical ejection. Considering 



these open questions. 



De Marcol (120091 ) reviews the possible shaping agents, remarking that 



binarity (whose definition is expanded to include 'normal' stars, brown dwarfs, and even 



planetary systems, as companions) is a necessary ingredient to explain n on-spherica . 



and possibly plays a fundamental role for spherical ones. In particular. 



De Marcol (120091 ) 



shap es 



considers common-envelope binary inte raction s as an original shaping agent in the formation 



of bipolar PNe, taking up the work of 



Sokerl ( 119971 ). Observations to probe binarity in 



central stars of PNe are non-trivial, but have been successful in a number of cases (e.g 



1979: 



Bond et al. 



Afsar fc Bond 



2009a 



1992 

liai 



Bond 



1995 



2000 



De Marco et al. 



2004 : 



Sorensen fc PoUaco 



20051) and rnore recently wit 



Frew fc Parker 



2007 



1 diffe rent techniques (JMiszalski et al. 



Miszalski et al. 



Bond 



2004 



star which is actually a close binary system. 



20081 . 



20111). A large fraction o f PNe have a central 



Miszalski et al 



( l2009bl ) note the fraction of 



close binary syst ems is ^ 17 ± 5 %, although there is evidence this number could be higher 



( iDe Marco 



20091). Currently, the number of close binary systems acting a s central stars of 



plane tary nebulae (CSPN) with orbital periods less than one day is ~40 (JMiszalski et al. 
201 ik 



Scenarios other than a close binary system have also been explored. Numerical 
simulations demonstrate it is possible to get an equatorial density enhancement in 
the superwind when stellar rotation and/or magnetic fields are included in the model 



-5- 



( JGarcia-Segura et al.lll999U2005l ). However, other stu dies have found that stellar differen tial 



rotation can be drained by the magnetic field itself ( ISoker 



2006 : 



Nordhaus et al. 



20071 ) 



reducing the field lifetime to ~ 100 yr, too short to influence the geometry of the gaseous 
ejecta. More work is needed to determine the extent to which stellar rotation and magnetic 
fields are important in the formation of single-star axisymmetrical PNe. 

To investigate the origin of PNe morphology and to characterize the important 
physical processes at work, theoretical models need to be supported with more high quality 
observational data. M orphokinem at ic studies are a fundamental technique available to 



approach the problem (JKwok 



20001 ). in which high resolution images and spectra are used to 
determine the real 3-D form of the complex gaseous structures. In this work we contribute 
to the study of PNe by modeling the morphokinematic structure of NGC 2818, an object 
whose highly detailed images have been published in the Hubble Heritage Project website.- 



1.2. NGC 2818 



Previous studies of NGC 2818 include the detailed work by iDufoun (Il984j ). 
who analyzed wide band imagery, photoelectric photometry, and optical and lUE 
spectroscopy. Dufour reports the following physical parameters: Te([0 III])= 14 500 ±500 K, 
Te([NlI])= 11 500 ± 800 K, iVe([S II])= 430 ± 250cm-3, and iVp([Cl III])= 1 300 ± 800 cm-^ 



Schildl fll995h obtained 



and a heliocentric systemic velocity of I^el = +18.5 ± 1.7kms~^. 

infrared images of NGC 2818, in which the central region was clearly seen. This region was 

more intense in the south, which was suggested to be related to shock excitation. 



Phillips &: Guest al (1l998[ l obtained narrow band Ha, [O III], [N II], and [S II] images 



from which some ratios were performed. In spite of the the low resolution of such images. 



^http://heritage. stsci.edu 
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the values of the ratios along the minor axis were plotted in diagnostic diagrams to 
investigate the nature of the emission. They did not directly observe clumps, but were able 
to deduce a clumpy medium formed by knots and filaments based on variations in electron 
density in the [S II] images. They also prop osed a possib le shock excitation for the southern 



region along the minor axis seen in H2 by 



Schildl ( 1l995l ). estimating planar shocks with 



veloci ties of V^ > 110 km s ^. Better resolution Ha and [O III] ir nages fromiGornv et al. 



(I1999I ) partially show the complex internal structure proposed by 



Phillips fc Cuestal f ll998h . 



A recent study in the Spitzer m id-infrared (MIR) bands (3.6, 4.5, 5.8, and 8/im) by 



Phillips Sz Ramos-Larios 



(I2OIOI ) shows gradients in the emission ratios. The observed 
gradients indicate an increasing importance of emission from photodissociation regions over 
emission from photoionization. Other contributions to the observed MIR emission include 
emission from polycyclic armatic hydrocarbons (PAHs) and effects of shocked regions. 

Finally, the study of NGC 2818 is of particular interest because it has been the topic 
of a series of articles in which its membership to a stellar open cluster is debated. This is 
an important issue given that distance is a parameter difficult to estimate for PNe, but 
relatively easy to estimate for stellar clu sters. In addition, there are only a few cases of PNe 



physically related to stellar clusters (e.g. 



Parker et al. 



20111 ) . Addressing this issue is outside 



the scope of this paper, but a detailed study to address this question is in preparation. 
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2. Observations and results 

2.1. Imaging 

A series of Hubble Space Telescope (HST) narrow band images, taken with the Wide 
Field Planetary Camera 2 (WFPC-2), were retrieved from the MAST Archive! to study 
the apparent morphology of NGC 2818 (proposal ID: 6119; PI: H. E. Bond; Date of 
observation: 1995 August 26). We show these images in Figured! presenting images with 
different contrasts for [N II]A6583 emission (filter F658N) in the two upper panels, Ha (left) 
and [O III]A5007 (right) emission lines (filters F656N and F502N) in the central panels, 
and [S II]A6731 (left) and He IIA4686 (right) emission lines (filters F673N and F469N) in 
the bottom panels. Three drizzled images were combined for each filter in order to remove 
cosmic rays and bad pixels according to standard procedures in HST reduction manuals. 

The appearance of NGC 2818 is slightly different when observed in different filters, 
as expected due to ionization stratification around the central star of a PN. In Fig. [H the 
basic structure appears to be composed of two bipolar lobes, separated in the center by a 
slightly pinched waist. The enhanced brightness in the cusps of the waist is similar to the 
limb-brightness of a denser squeezed ring-like structure. We will refer to this area around 
the central star, clearly seen as the unique emission region in He II, as 'the central region'. 

Both the bipolar lobes and the central region are visible in all panels, with the 
exception of He II, in which only the central region is visible. We note that the [O III] 
emission appears to be a smooth distribution of gas, filling the entire nebula. The length of 



^This study is based on observations made with the NASA/ESA HST, obtained from 
the data archive at the Space Telescope Science Institute (STScI). STScI is operated by 
the Association of Universities for Research in Astronomy, Inc. under NASA contract NAS 
5-26555. 



the major axis of the nebula projected on the plane of the sky is ^ 101 arcsec, whereas the 
minor axis is 31 arcsec, with the smallest separation of 14 arcsec at the cusps of the waist. 
The whole nebula can be contained in an ellipse of 126 x 62 arcsec. 

In addition to the main structure, a series of microstructures (cometary knots) are 
especially visible in the [N II] image, but can also be seen in the [S II] and Ha images 
(Fig. [H left panels). The tails of these cometary knots are primarily directed outwards from 
the geometrical center of the nebula. One of these microstructures, which is very prominent 
in the central region of the nebula, appears as an asymmetrical and broken 'big tail', almost 
aligned SE-NW. The upper-right panel of Fig. [T] shows a low contrast [N II] image in order 
to highlight the surroundings of the main structure. At the edges, several tails are also 
detected that can be traced outwards from the nebula. 

We made ratio maps using the HST images (Fig. [1]) to emphasize different morphological 
features in NGC 2818. Fig. [2] shows a mosaic of four image ratios. The upper-left panel 
shows the [N II] /Ha ratio in which the low-ionization regions (white) are better defined 



than in the original images. The image ratio [O III 
ionization/excitation and chemical abundances (see 



/Ha is sensitive to the variation of 



Corradi et a. 



19961) and can also be 



20071 ). Given the nature 



used to look for signatures of collimated outflows (iMedina et al. 

of the involved ions, the [S II] /Ha image ratio is sensitive to shocked regions whereas the 

He II/Ha image ratio shows the very high ionization regions. 

The location of the central star of NGC 2818 can be easily found by looking for a hot 
star in a blue-to-red image ratio. Fig. 12] shows an enlargement of the central region of the 
nebula in the [O III] /Ha image in which the central star appears as a central white dot. 
The location of the central star is in agreement with the extrapolation of the radial tails 
of many knots (not shown), as expected if they originate from the identified central star 
position. 
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2.2. Longslit high-dispersion spectroscopy 



Hi gh- dispersion optical 



(MES; 



Meaburn et al. 



spectra were obtained with the Manchester Echelle Spectrometer 



120031) and the 2.1m (f/7.5) telescope at the OAN-SPMfl observatory 
during 2009 February 4-6. A SlTe CCD with 1024 x 1024 pixels was used as detector. The 
slit length is 6'5 and the width was set to 150-yum (2"). A 2 x 2 binning was used, leading 
to a spatial scale of 0'.'6pixel~^ and a spectral scale of 0.1 Apixel"^. This spectrograph has 
no cross dispersion, consequently a AA = 90A bandwidth filter was used to isolate the S?*'^ 
order covering the spectral range which includes Ha and [N II]A6583 emission lines. The 
spectra were wavelength calibrated with a Th-Ar arc lamp to an accuracy of ±lkms^^. 
The FWHM of the arc lamp emission lines was measured to be ~ 12kms^^. Spectra were 
obtained by positioning the slit across five different regions of the nebula, as shown in 
Figure H The position angles (PAs) for these slits are +89° (slit 1), +123° (slit 2), and -1° 
(slits 3, 4, and 5). Exposure time for all the slits was 1800 s. Seeing was ~2 arcsec during 
observations. 

Fig. |5] shows the Position- Velocity (PV) maps for Ha (upper panels) and [N II] A6583 
(lower panels) emission lines. Each column is labeled with the corresponding slit number. 
We note that the He IIA6560 line is also observed in the PV maps for slits 1, 2, and 4, in 
agreement with the highest excitation location shown in the He II image (see lower-right 
panel in Fig. [T]). The relative position is measured with respect to the position of the central 
star (for slits 1, 2, and 4) and to the major axis defined by slit 1 (for slits 3 and 5). Radial 
velocity is relative to the systemic velocity I^el = +26 ±2 km s"""^ (Vlsr = +11 ±2kms^-'^). 
This value was deduced as the mean velocity of the two components of the splitting of 
the spectral lines, taken from the slit region crossing the central star (on slits 1, 2, and 



■^The Observatorio Astronomico Nacional at the Sierra de San Pedro Martir (OAN-SPM) 
is operated by the Instituto de Astronomia of the Universidad Nacional Autonoma de Mexico 
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4). This method only can be done if the sphtting of the spectral lines can be measured; 
that is, when we have high-resolution spectra. It has better accuracy than only taking 
the maximum of a broad unique emission line, as can be seen in low-resolution spectra. 
Observing without enough spectral resolution to see the splitting can result in a wrong 
estimate for the systemic velocity if, as in this case, the components of the splitting have 
different brightness. Without high resolution spectroscopy, the non-resolved splitting is 
seen as a broad emission line with its maximum being biased towards the wavelength of 
the brightest component. Previous determinations of the systemic velocity of N GC 2818 
have been been affected by low resolution spectra (e.g., 14ci = 18.5 ± 1.7kms~^ 



Dufour 



19841 ). In this paper, we present the first systemic velocity determination from the splitting 



method, utilizing high resolution spectroscopy. 

3. Discussion 

3.1. The morpho-kinematic structure of NGC 2818 

As described above (Sec. 12. ip . our inspection of the unprecedented high resolution 
HST images indicates that NGC 2818 is a morphologically rich PN with minor bubbles 
asymmetrical outflows, and cometary knots inside its main structure (Figs. [Hand [2]). 
In a first attempt to study the morphology and kinematics of NG C 2818, we have 



201l|),a 



constructed some simple models using the interactive software shape (iSteffen et al.l 
computational tool for the morpho-kinematic modeling and reconstruction of astrophysical 
objects. The internal library of shape includes basic structures such as spheres, cylinders, 
and tori, which can be modified by applying different parameters and/or functions 
("modifiers"). Each structure is composed of a grid in which a certain number of points are 
distributed. The user defines the number of points for each structure and the volume or 
surface distribution. Physical parameters, such as density and velocity, can be assigned to 
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the points either as an analytical functi on or interactive 
SHAPE to get the structure of PNe (e.g. 



Miszalski et al. 



y. Several recent 



2011 



Hsia et al 



papers have usee 



2010 



Jones et al. 



2010 



Contreras et al 



20101 ) ■ In our particular case, we obtained a final nebula model which 
fits the main structure of NGC 2818 for both imaging and spectroscopy (Fig. [6]). We find 
NGC 2818 can be modeled as a bipolar structure with major and minor axis length of 



1.8 X 0.3 pc (150 X 28arcsec at a distance of 2.5 kpc, following iVan de Steene fc Zijlstra 



19951 ). with an inclination angle of the major axis equal to 60° with respect to the line of 
sight (eastern lobe is blueshifted, western lobe is redshifted). The polar and equatorial 
velocities are Fpoi = 105kms~^ and Kq = 20kms~^, respectively, assuming homologous 
velocity expansion {V ex r). The shell has a maximum thickness of 0.1 pc (c^ 9arcsec) 
becoming thinner towards the poles. Using these velocities and dimensions, we have 
estimated a kinematical age of r^ ^ 8, 400 ± 3, 400 yr, which is simila r to the age of other 



PNe with cometary knots in their morphology (see lO'Dell et al. 



20021). 



The relatively high uncertainty in the kinematical age estimate merits a more detailed 
explanation. Distance estimates for PNe remain an important issue, given the very few cases 
in which we can get a 'direct' measurement (e.g. via stellar parallax). We have avoided 
distance estimates for NGC 2818 based on membership of this object to the stellar cluster 
mentioned in Sec. 11.21 and instead have used statistic al values . The re are some estimates 



for t 



le distance to this obiect r a nging from 1.79 kpc ((Phillips 



19951). 



20041 ) to 3.05 kpc flzhane 



Van de Steene fc Zijlstral (119951 ) give a value of 2.5 kpc, which is very close to mean 



value of the distance range. An uncertainty of ~40% is obtained utilizing the statistical 
estimates (which covers all the range of distance estimates). Thus, the propagation of this 
high uncertainty has effects in the same proportion over the physical size (1.8 ± 0.7 pc in the 
major axis) and the kinematical age (8,400 ± 3,400yr). Having better distance estimates 
will improve the accuracy on these two parameters, and is the focus of future work. 
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Because of their faintness, kinematics from the microstructures cannot be determined 
as they are not seen in our spectra, with the exception of the 'big tail' mentioned in Sec. 12.11 
A careful analysis of images in Figs. [H |21 and El as well as the PV map in Fig. E] (panel 2, 
lower) indicates that the 'big tail', as well as the other similar structures, have the same 
nature as the rest of the smaller cometary knots, being possibly photoevaporated, blown, 
and ionized by the action of the wind and radiation from the central star over the dense 
clumps that remain in its surroundings. In the 'big tail' a coUimated ejection is unlikely as 
the two NW pieces do not share the same kinematic behavior, as would be expected in the 
case of a collimated ejection (see Fig. \5\ panel 2, lower). If the 'big tail' is not taken into 
account, it is remarkable to note that the small cometary knots have no apparent preferred 
direction, whereas the denser and larger features are rela ted to the central reg ion. The 



19981 ). which has 



latter seems to be similar to the case of the Helix nebula ( iMeaburn et al. 

been proposed to have a knotty torus whose axis is seen at 37° with respect to the line of 

sight. 

Finally, the cometary knot antipodes in the upper-left side of the central region of NGC 
2818 appear to point towards a kind of 'local' bubble center instead of the central star. One 
possible interpretation is that the peculiar kinematics of several bubbles inside the nebula 
play an important role in the production of tails. We think that this idea deserves more 
attention in theoretical studies and is beyond the scope of the present work. 

3.2. The evolution of NGC 2818 

In this section we consider the line ratio images and what these ratios can tell us 
about the evolution of NGC 2818. First, we examine the image ratio [O III] /Ha. This 
ratio is sensitive to the variation of ionization/excitation and chemical abundances (see 



Corradi et al 



19961 ) and can also be used to look for signatures of collimated outflows 
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flMedinaet all 120071 ). According to Medina et all (120071 ) . PNe with [O III] /Ha ratios 
enhanced in an outer layer (or "skin"), are related to both, fast collimated outflows and 



multip 



Balick 



e-she lls, the latter may even have FLIERs (Fast Low-Ionization Emission Regions; 



19871). However, as seen in Fig. |2] (upper-right panel), NGC 2818 has a [O III] /Ha 



ratio diminished (Medina's type C), which could correspond to either young or evolved 
PNe, but not related to jets, shocks, or multiple shells. It is important to note that no 
correction for extinction was applied. In particular, this ratio image could be more affected 
by interstellar extinction (or internal extinction, including internal small-scale variatons) 
than the other ratios, given the relative large difference between the wavelengths of [O III] 
and Ha. 

The [S II] /Ha image ratio (Fig. [21 lower-left panel) is sensitive to the shocked regions 
and is enhanced in the most noticeable tails seen in [N II]. These could indicate either 
that these tails are shocked gas or that gas in the tails is being shocked by the stellar 
wind. Following the homologous velocity assumption, we expect no more than 28 kms~^ in 
stellar wind velocity at the distances in which the cometary knots are found, as they are 
preferentially located at a radius of 20arcsec around the central star (0.24 pc, assuming a 



distance of 2.5 kpc). However, these low 
[S II] /Ha ratio in a detectable way (see 



velocity shocks cou 



Phillips Sz Guzman 



d be enough to give rise to the 



1998 



and references therein). 



The same rat io also appea rs to exist in the ring zone. T his could be in agreement with 



the results of 



Schildl fll995h and 



Phillips &: Guest al (119981 ) who detect emission possibly 



related to shocked gas. It is also remarkable that based on the imaging in their study, 
Phillips Sz Guest al (119981 ) propose velocities of 110 km s^^ for the farthest regions of the 
nebula. Based on kinematics, we also estimate polar velocities around this value. 

The He Il/Ha image ratio shows the very high ionization regions. As expected, the 
highest ionized regions in NGG 2818 are enclosed inside the main structure of the nebula in 
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the regions surrounding the central star. This can also be seen in the the [O III] and He II 
images. In previous sections we defined the 'central region' based on the He II emission, 
whose appearance reminds one of the remnant of an equatorial structure. A possible 
interpretation of the data could be that NGC 2818 is a bipolar PN whose equatorial density 
enhancement can be traced by the cusps aligned with the central star and the brightness 
of the central region. The distorted internal structure and the presence of multiple bubbles 
and non-uniform microstructures can be seen as part of a pre-existent clumpy medium in 
the surroundings. Although a definite model for the formation of cometary knots has not 



been generally accepted, 



le.g- 



Capriotti fc Kendall 



lydro dynamical instabilities appear to adequately explain them 



20061 ). Such instabilities could be formed in the beginning of the 
proto-PN phase and appear as cometary knots once the bipolar outflow is free. Trying 
to understand the energetics involved in the whole process of formation and evolution of 
this nebula is an interesting way to address the study of the synergy between macro- and 
microstructures in the morphology of PNe. 

4. Conclusions 

We have conducted an observational study of NGC 2818 using archival high resolution 
space-based imaging and high dispersion Earth-based longslit spectroscopy. We derive 
a systemic velocity of Vhel = +26 ± 2kms^^. The computational software shape was 
used to construct a model to fit our data and improve our interpretation of the overall 
morphokinematic structure of NGC 2818. Our model is consistent with a bipolar nebula 
with a semi-major axis of 0.92 pc, possibly deformed by the stellar wind, a O.lTpc diameter 
central region that is a potential remnant of an equatorial enhancement, and a great 
number of cometary knots preferentially located inside a radius of 0.24 pc around the 
central star. The major axis of the main structure is oriented at z ~ 60° with respect 
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to the line-of-sight and at PA = +89° on the plane of the sky. Deprojected velocities 
correspond to an expansion of Vpoi = 105kms~^ and V^q = 20kms~^, which lead to an 
estimate of the kinematical age of r^ ~ 8, 400 ± 3, 400 yr. Cometary knots could be related 
to hydrodynamical instabilities formed before the bipolar outflow started. 

This project was supported by UNAM-PAPIIT grants IN111903 and IN109509. The 
author is grateful to the anonymous referee for his valuable comments and to his research 
group "Planetosos" for its support and fruitful discussions, in particular to Mr. Paco 
Beretta for calling my attention to this object. W. Steffen and L. Guiterrez provided SHAPE 
and HST advising, respectively. The author also acknowledges the OAN staff at San Pedro 
Martir, with special mention to Mr. Gustavo Melgoza-Kennedy (aka "Tiky") for support 
during observations. The author is deeply grateful to Dr. Ashley Zauderer for a kind and 
careful revision of the manuscript. Finally, the author would like to dedicate this paper in 
memory of Prof. Yolanda Gomez, who was an active researcher in the field of Planetary 
Nebulae and an indefatigable impeller of science education and outreach, and who recently 
passed away. 
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Fig. 1.— HST-WFPC2 narrow band images. Upper panels: [N II]A6583 (filter F658N) with 
two different contrasts; central panels: Ha (left; filter F656N) and [O III]A5007 (right; filter 
F502N); lower panels: [S II]A6731 (left; fiher F673N) and He IIA4686 (right; filter F469N). 
Spatial scale is indicated. Sky orientation (north up, east left) is the same in all subsequent 
figures. 
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Fig. 2. — HST-WFPC2 narrowband image ratios of NGC2818. White corresponds to large 
values of each ratio, indicated by the gray-scale bar. Upper left: [Nll]/Ha. Upper right: 
[Olll]/Ha. Lower left: [Sllj/Ha. Lower right: Hell/Ha. 
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Fig. 3. — Detailed view of the central region in the [O III] /Ha image ratio presented in 

Fig. 1. The central white dot corresponds to the central star of NGC 2818. 
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Fig. 4. — Slit positions used in the high dispersion spectroscopy superimposed on an [N II] 
HST image. PAs=+89° (1), +123° (2), and —1° (3, 4, and 5). The image was gaussian 
filtered to reproduce the spatial resolution observed in the spectra. 
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Fig. 5. — Position- Velocity maps of NGC2818 for slits positions 1 to 5. Upper panels 
correspond to Ha, and lower panels to [N ll] A6583. The slit orientation indicated in the 
upper panels is valid for the two maps in each column. The horizontal dashed lines represent 
the position of the central star for slits 1, 2, and 4, and the location where slits cross the 
major axis for slits 4 and 5. The horizontal dashed line corresponds to the systemic radial 
velocity (Vhel = +26 ± 2kms-i). 
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Fig. 6. — Comparison of observed data and the 3D model obtained from shape. Left: 
Observed [Nil] image (upper panel) and the shape image (lower panel). Right: Observed 
PV map of the [N h] A6583 line at slit position 1 (upper panel) and the corresponding shape 
synthetic PV map (lower panel). 



